Ginzburg-Landau theory of an RVB superconductor 
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We present a Ginzburg-Landau formulation of the bosonic resonating-valence-bond (RVB) theory 
of superconductivity. The superconducting order parameter is characterized by phase vortices that 
describe spinon excitations and the transition to the superconducting state occurs when such phase 
vortices (un)bind. We show that the boson RVB theory always leads to hc/2e flux quanta, and that 
the presence of a trapped spin-1/2 moment inside a vortex core gives rise to observable consequences 
for the low temperature field-dependent specific heat. We also show that the cores of magnetic 
fluxoids exhibit enhanced antiferromagnetic correlations. 
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Following the original proposal of Anderson |Q , there 
is, by now, a considerable body of literature devoted 
to the study of the resonating-valence-bond (RVB) the- 
ory. Much of this is based on the so called t — J model, 
which is the simplest possible model describing electrons 
(holes) moving in an antiferromagnetic background 
The RVB theories postulate that such a system is best 
described by collective spin and charge degrees of free- 
dom (usually called spinons and holons in the litera- 
ture), rather than a quasiparticle theory of interacting 
electrons. While there is no satisfactory proof that the 
quasiparticle theory fails in the t — J or related mod- 
els, there are some indications for such a failure, both 
from theoretical and numerical studies. Phenomenology 
of the cuprate superconductors based on these ideas has 
also been fairly successful. Motivated by these consider- 
ations, we present, in this paper, an effective theory of 
a superconducting RVB state. The theory leads to sev- 
eral interesting experimental consequences, and is very 
different from the conventional BCS theory of supercon- 
ductivity. Since the theory makes definite predictions, its 
veracity can be tested easily. 

The theory we present, is based on a bosonic descrip- 
tion of the t — J model. We choose this description since 
it accounts very well for the short-range antiferromag- 
netic (AF) correlations that play an important role in 
the cuprate superconductors. Further, in the limit of zero 
doping, the t — J model reduces to the Heisenberg model 
and the bosonic RVB theory in this limit (Schwinger- 
boson mean field theory) provides an excellent descrip- 
tion of the AF long-range ordered ground state, and the 
excited states ||, ||. In this limit, the bosonic RVB 
state can also be related to the variational wave function 
of Liang, Doucot and Anderson |J, as shown by Chen 
. Away from half filling, the theory describes bosonic 
charge (holon) and spin (spinon) degrees of freedom. In 
the bosonic theory of the RVB state , the electron oper- 
ator is expressed in terms of the bosonic holon and spinon 
operators and a topological (vortex) phase operator as 

c ia = h\b i(T e i§i ° . (1) 



ft satisfies usual fermionic anticommutation relations. 
The phase operator 0i CT is the most important ingredi- 
ent of the theory. It arises because the doped holes move 
in a spin background (not necessarily ordered) with AF 
correlations, and reflects the nonlocal effects of adding 
a hole to a doped Mott AF insulator. In what follows, 
we shall be concerned with the holon and spinon degrees 
of freedom, and therefore, it is natural to ask what role 
does the phase operator play in the description of these 
degrees of freedom. As shown in j7| , rewriting the t — J 
model in terms of the bosonic representation (fjl) leads 
to the emergence of two link fields, A s and A . The 
link field A s is coupled to the holon degrees of freedom 
and describes fictitious fluxoids bound to spinons satisfy- 
ing J2 c Afj = ±7r5^ Jec 0'n{ o ., for an arbitrary closed path 
c. Here, n\ denotes the spinon number operator. Simi- 
larly, the link field is coupled to the spinon degrees of 
freedom and describes fictitious fluxoids bound to holons 
satisfying J2 C A^ = ^Y^Uc n l' wnere nh 1S the holon 
number operator. Thus, we are led to the following phys- 
ical picture. The motion of holes in an AF background 
leads to nonlocal correlations between the charge and 
spin degrees of freedom, that are described by the link 
fields. The holons feel the presence of the spinons as vor- 
tices (quantized flux tubes) and vice versa (see Fig. 1). 
Inasmuch the holons and spinons perceive their mutual 
presence through the ir flux quanta, this theory can also 
be thought of as a tt flux theory, albeit one where the 7r 
fluxoids are bound to the constituent particles. 

The effective Hamiltonian is given by H e g = + H s , 
where and H s are the holon and spinon Hamiltonians 
respectively. Let us begin with the holon Hamiltonian, 
Hh- For convenience, we work with a continuum model 
and write 

Hh~7^— I d 2 v h\r) HV - A e - A s ) 2 h(r) , (2) 
2m/, J 

where rrih ~ (2ti l a 2 )~ 1 (th ~ t) is the effective mass of 
the holon, and A e , the vector potential of the external 
electromagnetic field. Note that the holon Hamiltonian 
(H) is coupled to the spinons through the term, A s (r). 
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This term is the continuum version of the link field A\ 
and is given by 



A*(r) = i 



r — r' 
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[n b T (r')-nJ(r')] . (3) 



The spinon Hamiltonian H s is given by 



J 

2 ^ 

<ij>a 



^ fa (<">■<■■■) hli,\ 



h.c. 



(4) 



where the spinon pairing (RVB) order parameter is de- 
fined by Af. = J^a < eT iaAii ^i^j-a >• Note again, 
the presence of the holon link variables in the definition 
of the spinon pairing (RVB) order parameter. For zero 
doping, the order parameter is identical to that of the 
Schwinger-boson mean field theory. As shown in |t]] , the 
Hamiltonian (Q) can be solved for a fixed hole concen- 
tration, within a self-consistent mean field theory. The 
bosonic spinon excitations are gapped, and the gap van- 
ishes as the hole doping decreases to zero. The short- 
range AF correlations are determined by the RVB order 
parameter Af-. 

When the holons (Bose) condense, the system should 
become superconducting, since the spinons are already 
paired. However, the presence of the phase field 0^ in 
(|I|) leads to interesting consequences. To see this, let us 
first write down the superconducting order parameter in 
terms of the decomposition (jl]). Assuming singlet pair- 
ing of electrons on nearest neighbor sites, we get (in the 
continuum limit), 



A^(r) 



A e i*'(r) 
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where A2 = f^A s [ip* h } 2 ( fa = ±1 for rj = x,y), 
■0 ; *(r) =< fv(r) > denotes the Bose condensate of the 
holons, and <J> s (r) arising from the phase field Oi CT in (|l|) 
is given by 



$ s (r) = / d 2 r' Im In [z - z'} (n\{r') - n\{r')) , (6) 

where z = x + iy. From (§) and (§), it is clear that 
<J> S describes phase vortices centered around the spinons: 
$ s — > $ s ± 2tt, if the coordinate r winds around a spinon 
continuously in space. Evidently, <& s is related to A s 
(V$ s = 2A S ) and both describe the vortex effect re- 
lated to the spinon excitations. From (^) and (pi), we see 
that a generalized Ginzburg-Landau (GL) equation can 
be written down for the holon condensate VVi(r) as 

aijj h + p\il> h \ 2 1> h + ^— HV - A e - A s ) 2 fl> h =0. (7) 
zmh 

The current operator can be constructed in the usual 
manner as 
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At T — and in the absence of an external magnetic 
field, all the spinons are (RVB) paired. Then, we see 
from (H) that A s is trivial and the phase of the super- 
conducting condensate is robust. But excited spinons can 
(and do) influence the holon condensate through the link 
field A s (r) in the above GL equation. In the following, 
we shall investigate such effects systematically. 

Superconductivity, confinement, and T c : From the def- 
inition of the order parameter (|^), we infer that the pres- 
ence of a holon condensate is not sufficient for supercon- 
ductivity to occur; i.e., phase coherence will not be real- 
ized in the system if the spinon vortices unbind, such that 
< e** ( r )e~ l * ( r ) > falls off exponentially This leads us 
to define the transition temperature, T c , as the temper- 
ature at which the spinon vortices bind, resulting in the 
vanishing of <i> s . To relate T c to the number of excited 
spinons, let us first consider the energy it costs to cre- 
ate an isolated spinon vortex. For a single spinon vortex 

X z x r 

centered around the origin, we have A s (r) = — — - — , 

2 r z 

for distances r >> a c , the size of the vortex core. From 

(||), we see that the presence of a spinon creates a su- 

percurrent in the holon condensate through A s [see also 
Ph 



(P) below]: J = — f^A s 

. mh . 
Substituting the above in > 



as illustrated by Fig. 2(a). 
3), we get 
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where L is the size of the sample. In general, the cost 
of creating more than one spinon vortex can be obtained 
by substituting the expression for A s (r), viz., (||) in (||). 
It is easily seen that the resulting expression describes 
a logarithmic attractive (repulsive) interaction between 
vortex-antivortex (vortex) excitations, like in the conven- 
tional Berezinskii-Kosterlitz-Thouless (BKT) transition. 

From the above, we conclude that a single S = 1/2 
bosonic spinon excitation is forbidden (owing to a log- 
arithmically diverging energy) and all excited spinons 
should be bound in vortex-antivortex pairs like in the 
low temperature phase of a BKT system. Consequently, 
$ s in (||) becomes trivial and phase coherence of (||) is re- 
alized. The ground state is therefore a superconductor in 
which 5=1/2 bosonic spinon excitations are confined by 
a logarithmic potential. However, such a confinement po- 
tential does not preclude the existence of fermionic quasi- 
particle excitations. Elsewhere, we showed that S = 1/2 
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nodal quasiparticles can be created from the condensate 
as a composite excitations. A calculation of the spec- 
tral function in the superconducting RVB state shows the 
presence of quasiparticle peaks below T c , consistent with 
results from angle-resolved photoemission spectroscopy 
M. Furthermore, an S = 1 spin excitation can be con- 
structed from a pair of bosonic 5 = 1/2 spinons, as shown 
in Fig. 2(b). The excitation energy E g , in this case, is 
finite and can be determined from the bosonic mean field 
theory [Q . This excitation leads to a sharp feature in the 
dynamical spin correlation function below T c . Note that 
the vorticities of the current vortices are not directly re- 
lated to the spin polarization directions, due to a gauge 
freedom to be discussed below. 

The transition to the normal state occurs when the 
phase coherence of the order parameter is destroyed; i.e., 
$ s in (||) is disordered owing to the emergence of free 
spinon vortices and the deconfinement of bosonic spinons 
marks the transition. It should be noted that the temper- 
ature at which the unbinding of spinon vortex-antivortcx 
excitations takes place can be substantially lower than 
the conventional BKT transition temperature (~ 1000 
K at optimal doping if A s is neglected ||). This is be- 
cause the cores of the spinon vortices begin to touch each 
other before the unbinding of vortices-antivortices driven 
by entropy happens. In this dense limit, the vortex- 
antivortex excitations can unbind because the energy cost 
is no longer logarithmically divergent, and T c is the tem- 
perature at which the average distance between spinons 
I ~ 2a c , where a c is the core radius of a spinon vortex, as 
illustrated in Fig. 3. The latter can be estimated by solv- 
ing the GL equation for an isolated vortex |l0fl . Here, we 
merely point out that the spinons have a characteristic 
length scale since their dynamics is governed by a uniform 
flux of Sir per plaquette. Recall that the spinons perceive 
the presence of holons as it fluxoids. The gauge field A^- 
in (|l|) represents an average flux of Sir (per plaquette) 
when the holons are condensed and the spinons acquire 
a characteristic length scale, a c ~ a/y/irS. Now, if n e s x 
be the number of excited spinons, then, I — 2a/ ^/Trn e s x . 
Therefore, we conclude that T c is determined by the 
condition n ex = 5. The RVB mean field theory of (^) 
shows that this condition is satisfied at a temperature 
T c ~ Eg /4 , Eg being the characteristic spinon energy 
0, as plotted in Fig. 3, where the experimental data 
from YBCO |ll[ are also shown for comparison (taking 
J = 100 meV). Such an estimate of T c based on the 
"core touching" mechanism agrees very well with the re- 
sults from a systematic renormalization group analysis of 
(§) by M. Shaw et al. @. 

Meissner effect and flux quantization: We now consider 
the situation when an external magnetic field is present 



Following the usual arguments for single valuedness of 
^(r), we get 



at T = 0. Writing ip h (r) = 
supercurrent given by (|8|) as 
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we express the 
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j>3{v)-dr = 2irn- jdr ■ (A e + A s ) 



(12) 



where the integral is over a closed loop and n, an integer. 
Now suppose that the integration is carried over a loop 
that is far away from the core of the vortex. Then, J = 
along the loop and we get 



27m - d> dr ■ A e - d> dr ■ A s = 



(13) 



When A s = 0, we see that the magnetic flux is quantized 
at 2wn in units of hc/e; i.e., the minimal flux quantum 
in this case is hc/e = 3>0j &s expected for a charge e 
Bose system. However, the presence of A s changes the 
quantization condition radically. Suppose there is one 
excited spinon trapped in the core of a magnetic fluxoid 
[Fig. 4(a)]. Then, from @, we obtain the minimal flux 
quantization condition, 



dr ■ A e = ±7T , 



(14) 



which is precisely the quantization condition of $o/2 m a 
superconductor with 2e pairing. As the holons do not dis- 
tinguish between internal (fictitious) and external (mag- 
netic) flux in (|2|) , they still perceive a total flux quantized 
at $o [see Fig. 4(a)], even though the true magnetic flux 
quantum is 3>o/2- 

Stability of the 2e flux quantum: When the magnetic 
flux quantum inside the core is <!>o, there is no spinon 
trapped inside the core. On the other hand, a bosonic 
spinon is trapped inside the core if the flux quantum is 
$0/2. Therefore, we have to estimate the energy dif- 
ference between these two cases to determine which of 
these is energetically favorable. The energy difference 
(per unit length) Ae between $o/2 and $ magnetic 
fluxes due to the magnetic field and supercurrent is given 

^— ^ J InK where A = (m^c 2 / Aire 2 5) , 

ft = X/a c , and a c ~ l/S 1 / 2 , as discussed earlier. Then, 

37T 

Ae = thShiK > —tS, which vanishes as S — > 0. So 

8 ~ 

far, we have ignored the cost of the vortex core. Since the 
$o/2 flux has a spinon trapped inside the core, we need to 
add the energy cost of an excited spinon to Ae estimated 
above. Since the energy of an excited spinon in the bulk 
(within the boson RVB theory @]) is E s = E g /2 ~ J 5, 
and since t > J, we expect the $0/2 flux quantum to be 
favorable for all doping concentrations p^ ] . Furthermore, 
as will be argued below, spinon (RVB) pairing inside the 
core is actually improved over the bulk; viz., the cost 
of creating a spinon inside the core is smaller than the 
energy of an excited spinon in the bulk. This further en- 
sures the stability of the <&o/2 flux quantum. This result 
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has to be contrasted with the slave boson RVB theories, 
where a $o flux quantum can be stabilized for small dop- 
ing O], In a recent paper |jl6), Wynn et al, carried 
out a careful search for $o fluxoids in the high temper- 
ature superconductors. Using scanning SQUID and Hall 
probe studies, they looked for <E> fluxoids in underdoped 
YBCO, but only observed <J> /2 fluxoids @. If spin 
charge separation indeed occurs in the high temperature 
superconductors, these experiments support our finding 
that only $o/2 fluxoids are favored energetically. 

Local moment inside a vortex core: As discussed above, 
each magnetic flux hc/2e is associated with an unpaired 
spin S — 1/2 trapped inside the vortex core (Fig. 4), 
where the holon density is suppressed. Note that the 
flux quantization condition (Q) is independent of the 
spin polarization. While <f c dr ■ A s = ±w does depend 
on the spinon polarization [c/.(0)], the resulting sign can 
always be absorbed by the phase of the holon condensate 
($ dr ■ V4>h = ±27r), without changing § dr ■ A e in (14). 
To paraphrase, each magnetic flux can trap either an 
S z = +1/2 or an S z = —1/2 spinon without changing 
the minimal quantization hc/2e. However, this two-fold 
degeneracy will be split by a Zeeman coupling to the 
magnetic field. This leads to a very unique consequence 
of the theory, discussed below. 

The Zeeman splitting between the S z = ±1/2 inside 
a vortex core can give rise to a Schottky contribution to 
the low temperature specific heat. From our analysis, we 
expect the number of magnetic moments [Fig. 4(b)] to 
be directly proportional to the applied magnetic field H, 
at least for small fields. The constant of proportionality 
can be estimated easily. Assuming a lattice constant of 
AA° for the Cu-0 planes, we find the number of induced 
magnetic moments given by 

^moment — TJ X H , 

where r\ = 8 x 10 -5 per Cu Tesla. Therefore, for a field 
of 6T, we estimate the number of induced magnetic mo- 
ments to be ~ 0.05% per Cu. Evidence for magnetic 
moments can be seen in the measurements of field depen- 
dent specific heat at low temperatures. However, most of 
these measurements are done on YBCO and the moments 
are invariably ascribed to impurities associated with the 
chains. Here, we offer another explanation based on our 
theory. The low temperature specific heat of the 123, 124 
and 247 phases of YBCO has been measured by Sanchez 
et al. ]l7| . They find that the observed specific heat in 
all these compounds can be described satisfactorily if one 
were to include a magnetic Schottky contribution whose 
amplitude (i.e., the number of magnetic moments) in- 
creases linearly with the field. For an applied field of 
6T, they estimate the number of induced magnetic mo- 
ments in all three phases of YBCO to be 0.07 - 0.08% 
per Cu atom. Similar results are observed in LSCO [Q, 
where the number of induced magnetic moments is also 



found to increase linearly for small fields (up to IT) be- 
fore saturating. Mason et al, [Q estimate the number 
of induced magnetic moments to be around 0.07% per 
Cu at 6T. These observations are in agreement with our 
estimates for the number of induced magnetic moments, 
as well as the field dependence. In another study of field 
dependent specific heat in YBCO, Emerson et al, |is| ] 
determined the number of induced magnetic moments to 
be 0.05 — 0.1% per Cu in the same range of magnetic 
fields. But they also find that the number of magnetic 
moments does not vary with the applied field. However, 
as the authors point out, their data is inherently different 
from the data of Sanchez et al. This discrepancy needs 
to be resolved by future measurements. Recently, Moler 
et al, 1 20 1 measured the field dependent specific heat of 
YBCO and also observed the number of magnetic mo- 
ments to be ~ 0.1% per Cu. 

Finally, we comment on the recent specific heat mea- 
surements in Zn doped YBCO pl| . These measurements 
show that the free magnetic moments (induced by Zn 
substitution) that are present in the normal state do not 
contribute to the specific heat below T c . There is no 
Schottky anomaly due to the magnetic moments for a 
wide range of Zn concentration in fields up to 8T. The 
authors explain this in terms of Kondo screening in the 
superconducting phase. However, the present work may 
provide a different explanation based on the confinement 
of bosonic spinons in the superconducting state. Ex- 
cept for the free moments trapped at vortex cores, no 
local moments related to S = 1/2 bosonic spinons are 
allowed in the superconducting bulk. This explains why 
the Schottky anomaly remains unchanged as in the Zn- 
free case up to 1% Zn concentration. 

Enhanced AF correlations: Thus far, we have seen that 
a spinon is trapped inside the core of a vortex. Clearly, 
the amplitude of the holon condensate ph vanishes at 
the trapped spinon position according to the GL equa- 
tion (Q), consistent with the single occupancy constraint. 
Since the spinon has a characteristic size a c , we expect 
the holon density to be reduced within such a core re- 
gion. This, in turn, enhances the RVB pairing inside the 
core. Consider the Hamiltonian (jij) for the spinons. At 
finite doping, the AF correlations are suppressed by the 
vortex effect of the holon link field A^.. Now, if the holon 
density is reduced inside the core region, then there is a 
reduction in the effect of A^ inside the core, and conse- 
quently, RVB pairing or equivalently, the short range AF 
correlations are enhanced compared to the bulk. In the 
extreme case of A^ — ► 0, when there are no holon vor- 
tices, the correlations inside the core will mimic those of 
the undoped compound. Thus, when a spinon is trapped 
near the vortex core, there is a concomitant enhance- 
ment of AF correlations inside the core region. This is a 
very important consequence of our theory and should be 
contrasted with the results from slave-boson RVB the- 
ories. In the latter, both the holon condensate and the 
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RVB pairing are suppressed inside the core region. Con- 
sequently, AF correlations are absent in the core region 
which is expected to be a normal Fermi liquid [0 , [Gl| . 
On the other hand, our considerations lead to the con- 
clusion that the core region is closer to the underdoped 
regime with enhanced AF correlations. 

Based on these arguments, we can obtain some sim- 
ple estimates. Recall that in the bulk, the characteristic 
length scale of a spinon, a c — a/VnS- Within the mag- 
netic vortex core, due to the reduction of S, we ex pect t he 
core size to increase, and the new core size, l c — \J8/S c a c , 
where S c denotes the average holon density within the 
core. To be consistent with the condition that one holon 
is expelled from the core region due to the trapping of an 
extra spinon, we demand p c x ttI 2 = p x nl 2 — 1, where 
Pc = 8 c /a 2 and p = S/a 2 . We then get S c = 5/2. Thus, 
the holon density is reduced by half within the core and 
l c = \[2a c . Since E g cc 5 in the bulk, we conclude that 
the characteristic spin excitation energy should be also 
approximately reduced by half, viz., E c ° Te ~ E g /2. Re- 
cently, Lake et at, reported the observation of AF corre- 
lations by neutron scattering inside the vortex core p2[ , 
where a field-induced sub-gap spin excitation was found 
in LSCO at an energy scale ~ 4.3 meV, approximately 
half of the characteristic spin energy scale in the bulk. 
The reduction of E c ° re due to the reduced holon den- 
sity may provide an explanation for such an observation. 
Tunneling spectroscopy of the vortex cores shows that 
the quasiparticle gap inside the core region is larger com- 
pared to the magnitude of the bulk superconducting gap 
p3[ . Again, this result does not contradict our conclusion 
that the core region is closer to the underdoped regime. 
These issues will be discussed in detail elsewhere p0[ . 

To conclude, we presented an effective theory of the su- 
perconducting state based on the bosonic representation 
of the t — J model. The theory leads to several interesting 
consequences. We first showed that the bosonic RVB the- 
ory leads to phase vortices in the superconducting order 
parameter. The phase vortices are excited spinons and 
T c is the temperature at which the cores of these vortices 
begin to overlap. We showed how hc/2e flux quantiza- 
tion leads to the trapping of a spinon inside the vortex 
core and argued that the cores exhibit enhanced anti- 
ferromagnetic correlations. Our estimates for the core 
energy show that the bosonic RVB theory does not al- 
low for hc/e flux quanta for any doping concentration. 
We also showed that the trapping of a spinon inside the 
vortex core leads to observable consequences for the low 
temperature specific heat. In a forthcoming publication, 
we shall present a quantitative analysis of the structure 
of an isolated vortex, based on the ideas outlined in this 
paper. We believe that the approach presented in this 
paper can bridge the gap between a microscopic model 
such as the t — J Hamiltonian and effective theories of 
quasiparticles in a d— wave superconductor coupled to 
fluctuating vortices p3. 
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Fig. 1. Holons and spinons perceive each other as 
carrying fictitious tt flux tubes. 

Fig. 2. (a) An S = 1/2 bosonic spinon carries a current 
vortex in the superconducting phase; (b) An S = 1 spin 
excitation is constructed from a pair of confined 5=1/2 
spinon vortices. 

Fig. 3. Phase coherence disappears at T c , as spinon 
vortices unbind due to core touching; T c is determined 
by the characteristic energy, E g , of the S — 1 excitation. 
The experimental data are from Ref. [ 111) . 

Fig. 4. Flux quantization occurs at hc/2e, with a 
bosonic 5=1/2 spinon trapped inside the core. 
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